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Abstract 

The small mass of the bottom quark, relative to its weak isospin part- 
ner, the top quark, makes the bottom an effective probe of new physics 
in Higgs and top sectors. We study the Higgs boson production associ- 
ated with bottom quarks, pp/pp — ► <j)bb — > bbbb, at the Fermilab Teva- 
tron and the CERN LHC. We find that strong and model-independent 
constraints on the size of the 4>-b-b coupling can be obtained for a wide 
range of Higgs boson masses. Their implications for the composite Higgs 
models with strong dynamics associated with the third family quarks 
(such as the top-condensate/topcolor models with naturally large bot- 
tom Yukawa couplings), and for the supersymmetric models with large 
tan (3, are analyzed. We conclude that the Tevatron and the LHC can 
put stringent bounds on these models, if the <pbb signal is not found. 
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I. INTRODUCTION 



A major task for all future high energy colliders is to determine the electroweak sym- 
metry breaking (EWSB) mechanism for generating the (W ,Z°) masses and the mech- 
anism for the fermion mass generation [3|. Whether the two mechanisms are correlated 
or not is an interesting and yet to be determined issue. Given the large top quark mass 
{m t = 175.6 ±5.5 GeV 0), as high as the EWSB scale, it has been speculated that the top 
quark may play a special role for the EWSB. One of such ideas is that some new strong 
dynamics may involve a composite Higgs sector to generate the EWSB and to provide a 
dynamical origin for the top quark mass generation (e.g., the top-condensate/top-color 
models ||). Another idea is realized in the supersymmetric theories in which the EWSB 
is driven radiatively by the large top quark Yukawa coupling || . 

In the minimal standard model (SM) there is only one Higgs doublet, which leaves 
a physical neutral scalar boson as the remnant of the spontaneous EWSB. The Yukawa 
couplings of the SM Higgs are determined from the relevant SM fermion masses divided by 
the vacuum expectation value (VEV), v ~ 246 GeV. Thus, aside from the coupling to the 
heavy top quark, all the other SM Yukawa couplings are highly suppressed, independent 
of the Higgs boson mass. For the top-condensate/topcolor type of models |§, with a 
composite Higgs sector, the new strong dynamics associated with the top sector plays a 
crucial role for generating the large top mass and (possibly) the W, Z boson masses. As 
to be discussed below, in this scenario, the interactions of the top-bottom sector with the 
composite Higgs bosons are different from that in the SM. Due to the infrared quasi-fixed- 
point structure || and the particular boundary conditions at the compositeness scale, the 
bottom Yukawa coupling to the relevant Higgs scalar is naturally as large as that of the 
top in such models. In the minimal supersymmetric extension of the SM (MSSM) ||, 
there are two Higgs doublets, whose mass spectrum includes two neutral scalars (h and 
H), a pseudoscalar (A) and a pair of charged scalars (H^. The MSSM Higgs sector 
contains two free parameters which are traditionally chosen as the ratio of the two Higgs 
VEV's (tan/? = v u /vd) and the pseudoscalar mass (m^). A distinct feature of the MSSM 
is that in the large tan /3 region the Higgs Yukawa couplings to all the down-type fermions 
are enhanced by either tan/3 or 1/ cos/3. Among the down-type fermions, the bottom 
quark has the largest mass and so the largest Yukawa coupling. Thus, it represents a 
likely place where new physics could reveal itself experimentally. This common feature, 
the large bottom Yukawa coupling relative to that of the SM, present in both types of the 
(conceptually quite distinct) models discussed above, serves as the theoretical motivation 
for our analysis. 
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Since the third family b quark, as the weak isospin partner of the top, can have large 
Yukawa coupling with the Higgs scalar(s) in both composite and supersymmetric models, 
we recently proposed to use the b quark as a probe of possible non-standard dynamics 
in Higgs and top sectors. In fact, because of the light b mass (~ 5 GeV) relative to 
that of the top (~ 175 GeV), the production of Higgs boson associated with b quarks 
(pp/pp — > <pbb — > bbbb) may be experimentally accessible at the Fermilab TevatronF] or 
the CERN Large Hadron Collider (LHC)Q, even though the large top mass could render 
associated Higgs production with top quarks (pp/pp — > <pti) infeasible. As we will show, 
this makes it possible for the Run II of the Tevatron and the LHC to confirm the various 
models in which the 6-quark Yukawa coupling is naturally enhanced relative to the SM 
prediction. However, if the <pbb signal is not found, the Tevatron and the LHC can 
put stringent constraints on the models with either a composite or a supersymmetric 
Higgs sector, in which the Yukawa coupling of the Higgs boson(s) and bottom quark is 
expected to be large. In H, this reaction was explored at the LHC and the Tevatron]] 
to probe the supersymmetry (SUSY) parameters of the MSSM. The conclusion was that, 
with efficient 6-tagging, useful information concerning the MSSM could be extracted from 
either the LHC or a high luminosity Tevatron. In this work, we expand upon earlier 
results J7[ in which it was concluded that even the much lower integrated luminosity of 
the Tevatron Run II can provide useful information through this reaction, provided an 
optimized search strategy is employed. We begin with a model-independent analysis in 
Sec. II, considering relevant backgrounds and determining an effective search strategy to 
extract a signal from the backgrounds. We then apply these results to constrain both 
composite and supersymmetric models in Sees. Ill and IV. In Sec. Ill, we first analyze 
the constraints on the two Higgs doublet extension of top-condensate model |J and the 
topcolor model ||10|| , and then remark upon the recent dynamical left-right model |11| 
(as a natural extension of the minimal top-condensate model [0) which also generically 
predict a large bottom Yukawa coupling. In Sec. IV, after deriving the exclusion contours 
on the m^-tan/3 plane of the MSSM, we further analyze its implication on the supergravity 



|]13fl and gauge-mediated |L4J models of soft SUSY breaking that naturally predict a large 
tan p. The comparison of our Tevatron Run II results with the LEP II bounds (from the 



1 A pp collider with y/s = 2 TeV. 

2 A pp collider with yfs = 14 TeV. 

3 In ||, a 1.8 TeV pp Tevatron was assumed with an integrated luminosity of 30 fb" 1 , and the 
squark mixings of the MSSM were neglected. 



2 



Zh and hA channels) is also presented, illustrating the complementarity of our analysis 
to other existing Higgs search strategies. Final conclusions are given in Sec. V. 



II. SIGNAL AND BACKGROUND 



We are interested in studying production of <pbb — > bbbb at the Run II of the Tevatron 
and the LHC. The signal events result from QCD production of a primary bb pair, with 
a Higgs boson (0) radiated from one of the bottom quark lines (c.f. Fig. [I]). The 
Higgs boson then decays into a secondary bb pair to form a bbbb final state. Because our 
detection strategy relies upon observing the primary b quarks in the final state (and thus 
demands that they have large transverse momentum), our calculation of the <j)bb signal 
rate from diagrams such as those shown in Fig. p] is expected to be reliable. This is in 
contrast to the inclusive rate of <fi production at a hadron collider, in which one does not 
require a final state topology with four distinct jets. In this case a calculation based upon 
Feynman diagrams such as those shown in Fig. |I| may not be reliable. It would be better 
to consider the Higgs boson production via bottom quark fusion, such as bb — > <f> and 



gb — > 4>b, with cares to avoid double counting its production rate |15|]. (This calculation 
would resum some large logarithms which are included in the definition of the bottom 
parton distribution function within the proton.) We have chosen to search in the four 
jet final topology because the QCD background for 3 jets is much larger than that for 4 
jets, and thus it would be more difficult to extract a 3 jet signal. Since the signal consists 
of four b (including b) jets, the dominant backgrounds at a hadron collider come from 
production of Zbb — > bbbb (c.f. Fig. 0), purely QCD production of bbbb (c.f. Fig. |) |16|] 
and bbjj, where j indicates a light quark or a gluon (c.f. Fig. |j) which can occasionally 
produce a 6-jet like signature in the detector. 

In order to derive model-independent bounds on the couplings of the scalar particles 
with the bottom quark, we consider K, the square-root of the enhancement factor for the 
production of (f)bb — > bbbb over the SM prediction. By definition, 

K y b 



(Vb)sM 

in which (yb)sM = V^rnb/v is the SM bottom Yukawa coupling and is the bottom 
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Yukawa coupling in the new physics model under the consider at ionF]. The decay branching 
ratio of <p to bb is model-dependent, and is not included in the calculations of this section. 
(Namely, the decay branching ratio Br(0 — > bb) is set to be one). We will properly take 
it into account for the specific models to be discussed in the following sections. 

We compute the signal and the backgrounds at the parton level, using leading order 
(LO) results from the MADGRAPH package fl7| for the signal and the backgrounds, 
including the sub-processes initiated by qq and gg (and in the case of bbjj, qg and qg). 
While the complete next-to-leading order (NLO) calculations are not currently available 
for the signal or background cross sections, we draw upon existing results for high p? 
bb production at hadron colliders |18[ and thus estimate the NLO effects by including a 
fc-factor of 2 for all of the signal and background rates. (NLO effects to the pp — > (pti 
cross section in the limit s 3> m\ 3> were explored in It was found that QCD 
fc-factor is on the order of 1.5, but this limit is not expected to provide a very good 
estimate to the <pbb rate at the Tevatron since the corresponding condition ml ^> m| 
no longer holds. Furthermore, because the mass of the bottom quark is much less than 
that of the top quark, we expect that the fc-factor for the <pbb production rate to be 
larger than that for the <frti rate.) In the end of this section, we will also estimate the 
uncertainty in the calculation of the signal and background cross sections based upon the 
above prescription. We use the CTEQ4L [f2(]] parton distribution functions (PDFs) and 
set the factorization scale, /io, to the average of the transverse masses of the primary b 
quarks, and the boson (0 or Z) transverse massQ for the 4>bb and Zbb processes, and use a 
factorization scale of //o = V§, where § is the square of the partonic center of mass energy, 
for the bbbb and bbjj background processes. It is expected that a large part of the total 
QCD bbbb and bbjj rates at the Tevatron or LHC energies will come from fragmentation 
effects, which we have neglected in our matrix element calculation. However, due to the 
strong pt and isolation cuts which we will impose for improving the signal-to-background 
ratio (explained below), we expect that these effects will be suppressed, and thus will 
only have a small effect on our results. Similarly, we expect that after imposing the 
necessary kinematic cuts, the signal and the background rates are less sensitive to the 
above choice of the factorization scale. In this section, unless otherwise noted, we will 



4 For simplicity, we ignore running effects in (y&)sM in this section, treating the bottom quark 
mass as 5 GeV at all scales. We will comment on the running effects of the Yukawa coupling in 
the context of specific models of new physics in Sees. Ill and IV. 

■ The transverse mass of particle , is gtven by „.<!> S -f^W . 
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restrict our discussion of numerical results to a signal rate corresponding to a scalar mass 
of 7«0 = 100 GeV, and an enhancement factor of K = m t /mb ~ 40. We will consider the 
experimental limits which may be placed on K as a function below. 

In order to simulate the detector acceptance, we require the px of all four of the final 
state jets to be pt > 15 GeV, and that they lie in the central region of the detector, with 
rapidity \r)\ < 2. We also demand that the jets are resolvable as separate objects, requiring 



a cone separation of AR > 0.4, where AR = y A(p 2 + Arj 2 . (A(p is the separation in the 
azimuthal angles.) In the second column of Table | we present the number of events in the 
signal and background processes at the Tevatron Run II which satisfy these acceptance 
cuts, assuming 2 fb _1 of integrated luminosity. As can be seen, the large background 
makes it difficult to observe a signal in the absence of a carefully tuned search strategy to 
enhance the signal-to-background ratio. In presenting these numbers, we have assumed 
that it will be possible to trigger on events containing high pt jets (and thus retain all of 
the signal and background events). This capability is essential for our analysis. 

The typical topology of the bottom quarks in the signal events is a "lop-sided" struc- 
ture in which one of the bottom quarks from the Higgs decay has a rather high p? of 
about 7770/2, whereas the other three are typically much softer. Thus, the signal events 
typically have one bottom quark which is much more energetic than the other three. On 
the other hand, the QCD bbbb (or bbjj) background is typically much more symmetrical, 
with pairs of bottom quarks (or fake fe's) with comparable pt- In order to exploit this, we 
order the b quarks by their transverse momentum, 

p?>p?>p?>p¥\ (2) 



and require that the bottom quark with highest transverse momentum have p T > 50 
GeV, and that p { ? > 30 GeV and p { ^ A) > 20 GeV.Q In Fig. |, we show the differential 
cross sections with respect to Pj? after the acceptance cuts and with respect to p^p after 
the cut for the signal and the QCD bbbb background at the Tevatron Run II. These 
plots illustrate the advantage in isolating the signal from the background provided by the 
asymmetric cuts on the final state b quarks outlined above. In the third column of Table | 
we show the effect of these cuts on the signal and backgrounds. As can be seen, these cuts 
reduce the signal by about 60%, while drastically reducing the QCD bbbb background by 



6 Here, we have corrected the values given in Ref. Q due to a Fortran error in evaluating the 
signal distribution. We thank S. Mrenna for cooperation in checking the Monte Carlo simulation 
and for his help in detecting this numerical error. 
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about 90%. 



Since the px spectrum of the leading jets is determined by the mass of the scalar boson 
produced, the leading p T cuts can be optimized to search for a particular m^. The results 
for several values of are presented in Table [T| As expected from the discussion above, 
the optimal cut on is close to whereas the optimal cut on is somewhat lower 

(generally closer to m^/3). We adopt these optimized pt cuts for each mass considered, 
when estimating the search reach of the Tevatron or LHC. 

Another effective method for reducing the QCD background is to tighten the isolation 
cut on the bottom quarks. In the QCD bbbb background, one of the bb pairs is preferentially 
produced from gluon splitting. Because of the collinear enhancement, the invariant mass 
of this bb pair tends to be small, and the AR separation of these two fe's prefers to be 
as small as possible. On the contrary, in the signal events, the invariant mass of the bb 
pair from the 0-decay is on the order of m^, and the AR separation is large because the 
angular distribution of b in the rest frame of the scalar <p is flat. Thus, by increasing 
the cut on AR to AR > 0.9 we can improve the significance of the signal. As shown in 
column four of Table |, this cut further decreases the signal by about 30%, and the QCD 
bbbb background by about 65%. In the end, their event rates are about the same. 

One can further improve the significance of the signal by attempting to reconstruct 
the mass of the scalar resonance. This can be difficult in principle, because one does not 
know a priori what this mass is, or which bottom quarks resulted from the (f) decay in 
a given event. It may be possible to locate the peak in the invariant mass distribution 
of the secondary b quarks resulting from the decay, though with limited statistics and 
a poor mass resolution this may prove impractical. However, one can also scan through 
a set of masses, and provide 95% C.L. limits on the presence of a Higgs boson (with a 
given enhancement to the cross section, K) in the bbbb data sample for each value of m$ 
in the set. In order to do this, we assume a Higgs mass, and find the pair of b quarks 
with invariant mass which best reconstructs this assumed mass. We reject the event if 
this "best reconstructed" mass is more than 2Am^ away from our assumed mass, where 
2Am ( j > is the maximum of either twice the natural width of the scalar under study (1^) 
or the twice experimental mass resolution^ As shown in the fifth column of Table |, this 



We estimate the experimental mass resolution for an object of mass to be Am^ = 



models of new physics considered in Sees. Ill and IV is usually smaller than this experimental 
mass resolution. 




Under this assumption, the natural width of the bosons in the specific 
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cut has virtually no effect on the signal or Zbb background (for a 100 GeV Higgs) while 
removing about another 10% of the bbbb background. 

As will be discussed below, the natural width of the Higgs bosons in both the MSSM 
and the models of strong EWSB that we wish to probe in this paper are generally much 
smaller than our estimated experimental mass resolution, and thus one might think that 
an improved experimental mass resolution could considerably improve the limits one may 
place on a scalar particle with a strong b interaction. However, the models in which we 
are interested generally have one or more nearly mass-degenerate bosons with similarly 
enhanced bottom Yukawa couplings. If the extra scalars are much closer in mass than the 
experimental mass resolution (and the natural width of the bosons), the signal can thus 
include separate signals from more than one of them. Thus there is potentially a trade- 
off in the AM cut (cf. Table I) between reduction of the background and acceptance 
of the signal from more than one scalar resonance. In order to estimate the potential 
improvement for discovering a single Higgs boson, we have examined the effect on the 
significance one obtains if the cut on the invariant mass which best reconstructs 
is reduced to Am^ as opposed to 2Am^ as was considered above. We find that this 
improved mass resolution further reduces the QCD bbbb background by about another 
40%. Assuming four b tags (as discussed below), this improved mass resolution increases 
the significance of the signal from about 12.2 to 14.6, which will improve the model- 
independent lower bound on K by about 10%. Thus, an improved mass resolution would 
most likely be helpful in this analysis. 

Another method to further suppress background rate is to observe that in the back- 
ground events, the b quarks whose invariant mass best reconstructs come from the 
same gluon. This is because, after imposing all the kinematical cuts discussed above (cf. 
Table I), the matrix elements are dominated by Feynman diagrams in which one very 
far off-shell gluon decays into a bb pair, as opposed to interference of many production 
diagrams, which dominates the lower invariant mass region. Thus, for greater than 
about 100 GeV, the background event produces b quarks with the characteristic angular 
distribution of a vector decaying into fermions, 1 + cos 2 9, in the rest frame of the bb 
system. This is distinct from the signal distribution, which comes from a scalar decay, 
and is flat in cos 9. Thus, for masses above 100 GeV, we further require | cos9\ < 0.7 after 
boosting back to the rest frame of the bb pair which we have identified as coming from 
the scalar boson 0. 

In order to deal with the large QCD bbjj background, it is important to be able to 
distinguish jets initiated by b quarks from those resulting from light quarks or gluons. 
We estimate the probability to successfully identify a b quark passing the acceptance cuts 
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outlined above to be 60%, with a probability of 0.5% to misidentify a jet coming from a 
light quark or gluon as a b jet ||21|| . In Table [TTT] we show the resulting number of signal 
and background events passing our optimized cuts at the Tevatron, assuming 2 fb _1 of 
integrated luminosity, after demanding that two or more, three or more, or four 6-tags 
be present in the events, and the resulting significance of the signal (computed as the 
number of signal events divided by the square root of the number of background events). 
We find that requiring 3 or more 6-tags results in about the same significance of 12.2a as 
requiring 4 6-tags. However, we see that for the chosen parameters (m^ = 100 GeV and 
K = rrit/miy ~ 40), even with only 2 or more 6-tags, one arrives at a significance of about 
3a, and thus has some ability to probe a limited region of parameters. From the large 
significance, we see that the Tevatron may be used to place strong constraints on Higgs 
particles with enhanced bottom quark Yukawa couplings, and that the ability to tag 3 
or more of the bottom quarks present in the signal can probe a larger class of models 
(or parameter space of the models) as compared to what is possible if only 2 or more of 
the bottom quarks are tagged. In the analysis below, to allow for the possibility that the 
bbjj background may be somewhat larger than our estimates, we require 4 6-tags, though 
as we have demonstrated above, we do not expect a large change in the results if 3 or 4 
6-tags were required instead. 

This analysis can be repeated for any value of m^, using the corresponding px cuts 
shown in Table [TT]. It is interesting to note that the signal composition in terms of the gg 
or qq initial state depends on the collider type and the mass of the produced boson, which 
controls the type of PDF and the typical region of x ~ vr?±jS at which it is evaluated. 
At the Tevatron, for = 100 GeV, the signal is 99% gg initial state before cuts, and 
87% after cuts, while for = 200 GeV, it is 99% gg initial state before cuts, and 85% 
after cuts. Thus, at the Tevatron, one ignores about 15% of the signal if one relies on 
a calculation employing only the gg initial state. At the LHC, for = 100, the signal 
is very close to 100% gg initial state before cuts and 99% after cuts, and for = 500 
GeV, it is 99% gg initial state before cuts, and 99% after cuts. This indicates that at 
the LHC, very accurate results are possible from a calculation considering only the gg 
initital state. The resulting numbers of signal and (total) background events after cuts 
for various boson masses are shown in Table [TV] . 

From these results, one may derive the minimum value of K, K m i n , for a scalar boson 
with mass to be discovered at the Tevatron or the LHC via the production mode 
bbcj){— > bb). Similarly, if signal is not found, one can exclude models which predict the 
enhancement factor K to be larger than K min . To give a model- independent result, 
we assume that the width of the is much less than the estimated experimental mass 
resolution defined above, which is the case for the models studied in this paper. We 
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determine K min by noting that in the presence of a Higgs boson with enhanced bottom 
Yukawa couplings, the number of expected signal events passing our selection criterion 
is given by N$ = K 2 Ng SM ^ , where Ng SM ^ is the number of signal events expected for a 
scalar of mass with SM coupling to the b quark [assuming Br(0 — > bb) = 1], whereas 
the number of background events expected to pass our cuts, N B , is independent of K. 
Thus, requiring that no 95% C.L. deviation is observed in the bbbb data sample (and 
assuming Gaussian statistics) determines 



K ■ 

J1 mm 



1.96 VN- B 



(SM) 



where 1.96a is the 95% C.L. in Gaussian statistics. In Fig. [| we show the resulting 
95% C.L. limits one may impose on K min as a function of from the Tevatron with 2, 
10, and 30 fb _1 and from the LHC with 100 fb _1 , as well as the discovery reach of the 
LHC at the 5a level. Our conclusions concerning the LHC's ability to probe a Higgs 
boson with an enhanced b Yukawa coupling are very similar to those drawn in M , but are 
considerably more optimistic than those in where the conclusion was that the bb jj 
background is considerably larger than our estimate (though there are elements of the 
search strategy which differ between those of p2| and ours as well, and their simulation 



of the ATLAS detector is certainly more sophisticated). In |22| the backgrounds were 



simulated using PYTHIA [23| to generate two to two hard scatterings and then generating 
the additional jets from a parton showering algorithm. As noted above, in the light of the 
strong (ordering of) pr and isolation cuts applied to select the signal events, we feel that 
a genuine four body matrix element calculation such as was used in our analysis provides 
a more reliable estimate of this background. 

We have examined the scale and PDF dependence of our calculation for the signal 
and background rates at the Tevatron, and find that in varying the scale between one 
half and twice its default choice (defined above), \i = Hq/2 and \i = 2/i , the (pbb signal 
and Zbb background rates both vary from the result at fi = //q by about 30%, while the 
bbbb and bbjj backgrounds vary by about 45%. This strong scale dependence is indicative 
of the possibility of large higher order corrections to the leading order rate. Thus, in 
order to better understand the true signal and background rates, it would be useful to 
pursue these calculations to NLO. We have also compared the difference in the results 
from the MRRS(Rl) PDF JZJ and the CTEQ4L PDF, and find a variation of about 10% 



in the resulting signal and background rates. Since these separate sources of uncertainty 
(from PDF and scale dependence) are non-Gaussianly distributed, there is no way to 
rigorously combine them. Thus, we conservatively choose to add them linearly, finding a 
total uncertainty of about 40% in the signal rate (Ng SM ^), and 50% in the background 
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rate (Nb). From the derivation of K min above, we see that these uncertainties in signal 
and background rate (which we assume to be uncorrelated) combine to give a fractional 
uncertainty in K min , 



5K ■ 
K ■ 

1 v mm 



SN ( S SM) \ fSX, 



\\2N { S SM) ) \4N 



B , 



(4) 



where SNg SM ^ and SNb are the absolute uncertainties in N$ SM ^ and Nb, respectively. 
From this result, we see that in terms of a more precise theoretical determination of 
K m i n , one gains much more from a better understanding of the signal rate than a better 
determination of the backgrounds. Applying our estimation of the uncertainty from PDF 
and scale dependence to Eq. we find an over-all theoretical uncertainty in K min of 
about 25%. 



III. CONSTRAINTS AND IMPLICATIONS ON DYNAMICAL MODELS WITH 
STRONGLY COUPLED ELECTROWEAK SECTOR 

The observed large top mass, of the order of the electroweak scale, singles out top 
quark from all the other light fermions. This makes the top-condensate/topcolor scenario 
particularly attractive In this section, we analyze the strongly interacting scenario of 
the EWSB with a composite Higgs sector. We consider top-condensate/topcolor type of 
models [p HT2]j25| -p8| in which new strong dynamics associated with the top quark sector 



plays a crucial role in the generation of the top quark and the W, Z boson masses. As we 
have emphasized, since the bottom quark is the weak isospin partner of the top quark, 
its interaction to the Higgs sector can be closely related to that of the top quark. In 
the top-condensate/topcolor scenario to be analyzed below, the b quark Yukawa coupling 
(to the relevant scalar) is naturally large, of the same order as the top Yukawa coupling 
[~ 0(1)], due to the quasi-infrared fixed point structure f| and the proper boundary 
conditions at the compositeness scale. This can give distinct experimental signatures at 
the Tevatron and the LHC In the following, we shall analyze two specific models in this 
scenario, and derive the constraints (or discovery reach) expected at the Tevatron and the 
LHC based upon the model-independent results in Sec. II. Finally, we shall analyze the 
dynamical left-right symmetric extension [[II] of the minimal top-condensate model [II 
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A. In Two Higgs Doublet Extension of the Minimal Top-Condensate Model 



Since the minimal top-condensate model (with three families) [T2[] was ruled out due to 
predicting a too large top mass (~ 220 — 250 GeV) to reconcile with experimental data, we 
consider the minimal two Higgs doublet extension (2HDE) of the top-condensate model 
proposed in Ref. @, which predicts a smaller value of m;. Though the simplest 2HDE of 
the top-condensate model may not provide enough reduction of the m t value to match the 
Tevatron measurement, it is not out of question to incorporate some further improvements 
H, e.g. including the recently proposed seesaw-type top-condensation PBfl, to achieve a 



realistic m t . In the supersymmetrized version of the top-condensate model P9 |, it is 
possible to derive the correct top mass while keeping the similar boundary conditions 
and the quasi-infrared fixed point structure which ensures the large b quark (and also 
r lepton) Yukawa couplings. In this subsection we analyze the simplest 2HDE of the 
top-condensate model constructed in Ref. 0, for illustration. 

The starting point of the model is to consider the SM without an elementary Higgs 
boson, but with Nambu-Jona-Lasinio (NJL) type of four- Fermi interactions [|3(J generated 
at the cut-off scale A, where the new physics takes place. For the third generation quarks, 
the SU(2)l (g> U(l)y invariant 4- Fermi couplings can be written as [|§: 



Uf = G t (* L t R ) (t R * L ) + G b (v L b R ) (b R * L ) + G tb 



h.c 



(5) 



where the summation over color indices is implied in the round parentheses and ^ = 
(—b c ,t c ) T . Then, just below the scale A, two composite Higgs doublets, $t and can 
be introduced via the auxiliary field method [[H]] with the interaction Lagrangian 



SF 



(y L $ t t R + h.c. 



,b R + h.c.) - i4 b + h.c. 



(6) 



To diagonize the $ t and $b mass terms, one needs to introduce the mixing angle a defined 
by: a = 2^ b /(fif — //f) . In Eqs. (||) and (f|), the mixing term proportional to G t b or 
is important to break the Peccei-Quinn U(l) symmetry and to generate a nonzero mass 
for the pseudoscalar. The low energy Lagrangian at the scale /i(< A) can be deduced 
from (|6|) via the renormalization group (RG) evolution JJIJ which defines the effective low 
energy couplings. Thus, at the scale fi , 

4f = Zlvt (^Mr + h.c.) + Zl b y b (* L <f> b b R + h.c.) + 

Z* t {D^ t )\D^ t ) + Z^ b (D^ b y(D^ b ) + V(Zl$ t , Zl b $ h ) , 
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where V is the renormalized Higgs potential for $ t and As fi approaches to A, (|7]) 
should match with the bare Lagrangian (|6]) to result in the proper boundary conditions to 
be used in the renormalization group analysis pi. It turns out that the compositeness of 
both $ 4 and can be achieved only for the boundary condition yt(A.) = y&(A) = y ^> 1 
II. Hence, 

Vt(fj) ~ Vb{^) , (for any [M < A) . (8) 



This is due to the fact that y t and yb satisfy similar RG equations except for a small 
difference (in the g\ term) originating from the different hyper-charges of the t and b 
quarks p2L§: 



1 



dybilA = 

tiln/x 167T 2 



16vr 2 
1 



| + y 6 V) + -yfM - 3 (jv 6 - ir) 53 2 (/^) - \al^) - ^<7 2 (m) 



(9) 



where N c = 3 for the QCD theory, and fi > m t ~ MHi ggs . Because of the infrared 
quasi-fixed point structure |p|,p2| of the two Higgs doublet model, the relation (^) holds 
well as long as y&(A),y t (A) > 1. This is true even for the case where yb is chosen to be 
significantly lower than y t at the compositeness scale A. This running behavior is shown 
in Fig. [j], which confirms the large value of yb at the weak scale. We have also examined 
possible threshold effects due to different values of the Higgs masses and found the above 
conclusion unchanged. 

The two composite Higgs doublets develop VEVs from the condensation of <it> 
and <bb> (determined by the gap equations), so that (v t ,0) T / \/2 and 

<$b>= (0,ffe) T / y/2 . Since the masses of the t and b quarks are given by m tj f,(/i) = 
yt,b(f-i>)vt,b(ti>) / V% ) an d the Yukawa couplings yb and y t are about the same, of the 0(1) 
at the weak scale [cf. (§) and Fig. 0], this model naturally predicts a large tan/?: 

v t {m t ) m t {m t ) 



tan (3 



v b (m t ) m b {m t ) 



55 > 1 



(10) 



Here, mb(t){ m t) is the running bottom (top) mass at the scale m t . The running values of 
ra&(/x) and m t (/x) are derived |H| from the measured physical pole masses ml° l ~ 5 GeV 
and m\ oX ~ 175 GeV, and are dominated by the QCD evolution at scales ~ 0(M H ig gs ). 



At the one- loop level, the relation between the pole quark mass m^° l and the MS QCD 
running mass at the scale /i 



mJ o1 is: 
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3tt 



'ID 



When running upward to any scale /i, 



m q (fi) = m q (m- 



pol\ 



c[a s (/i)/vr] 



a s (m q ml )/n 



(12) 



where c[x] = (23x/6) 12,/23 [1 + 1.175a;] for m£ o1 < \i < mf o1 , and c(x) = 
(7a;/2) 4 / 7 [1 + 1.398a;] for /i > mj? Dl |[35|| . Numerically, m t (mf o1 ) ~ 166 GeV and 



m h \m\ 



pol\ 



3 GeV. 



Analyzing the mass spectrum of the Higgs sector, we find that the lightest scalar 
particle with the large Yukawa coupling to the bottom quark is the pseudoscalar P (= 
y/2 [sin /3Im$g + cos /3Im$ t °] ~ y/2Im$° b ) with a mass 



V2 



A 4 | 2 



tan(vr - 2/5) 
tan 2a 



- 1 



(13) 



Since v = Jv% + t> 2 ~ 246 GeV, M P can be as low as 0(m z ), depending on the Higgs 
mixing angle a and the Higgs self-couplings A4 [f|. For instance, for A = 10 15 GeV, the 
mass Mp is less than about 233 GeV, and the decay branching ratio of P to bb is about 
one. Hence, this model predicts a light pseudoscalar that couples to the bottom quark 
strongly through Yukawa interaction. 

To discover or exclude this model at the Tevatron and the LHC via measuring the 
production rate of four b jets, we need to make use of the model-independent results 
obtained in Sec. II, Namely, we need to compare the predicted bottom Yukawa coupling 
[Vb = Vbil 1 ) sin 0\ of the pseudoscalar P with the model-independent bound (= K m i n y^) 
derived in Sec. II, where the reference value is arbitrarily chosen to be = y/2ml° l /v 
with m%° 1 ~ 5 GeV and v ~ 246 GeV. This is equivalent to comparing y^ /yf^ 1 with 
K min , where y b is the running Yukawa coupling at the scale Mp. Note that at the weak 
scale the running effects of the VEVs [mainly due to the electroweak corrections] are 
negligible and for fi ~ M P , m b (fi) = yb{fi)v b = y^ M (fi)v, so that yf (^) = sin /%(//) = 
tan/3 \/2mi,(fx) / v . We can thus derive a minimal tan/3 value for a given Mp to discover 
such a model at hadron colliders. The result of 95% C.L. exclusion contours is given in 
Fig. 8 for various colliders. As shown in Fig. |8|, the 2 fb _1 Tevatron Run II data can 
exclude models with Mp up to ~ 190 GeV, if no signal is found. The entire mass range of 
P predicted in this model (less than 233 GeV for A = 10 15 GeV) can already be explored 
at the Tevatron Run II with a 10 fb _1 luminosity. 



13 



B. In Topcolor Assisted Technicolor Model 



The minimal top-condensate model and its two-Higgs-doublet extension || re- 
quire fine-tuning the four-Fermi coupling(s) at the scale A to be very close to the 
critical value, but do not address the dynamical origin of the effective coupling(s) at 
the energy scale above A. The topcolor assisted technicolor models (TCATC) [fTOf] 
were proposed to overcome such difficulties. These models postulate a gauge structure 
g = SU(3)i ® SU(3) 2 ® £7(l)i ® U{1) 2 ® SU(2) W at the scales above A = 0(1 TeV) . 
The third family fermions couple to S77(3)i ® ?7(l)x gauge sector with the same quan- 
tum numbers as those under the SM QCD and U(l)y interactions, while the first two 
family fermions couple to SU(3)2 <§) f (1)2 m a similar way. At the scale A = 0(1 TeV), 
SU(3)i ® U(l)x is strong but not confining, and Q is spontaneously broken down to 
Gsm = SU(3) C <S> U(l)y <S> SU{2) w due to an unspecified mechanism which may or may 
not be related to the EWSB. In consequence, massive gauge bosons of the color octet B a 
(colorons) and the singlet Z' are generated. Below this breaking scale A = min(M^, Mz>), 
4-Fermi interactions are generated as follows: 



C-af — — 



4it 
A2 



2ki 



9N r 



+ 



9N r 



(14) 



After the fermions condensate, an effective Lagrangian with two composite Higgs doublets 
(<&t and $&) can be introduced as: 



C 



SF 



y t ^ L ^ t t R + y b V L $ b b R + h.c 



- A 



(15) 



with 



y t = J47t(k + 2/^/9^) 



y b = ^Ati(k- k x /9N c ) 



(16) 



Here k and K\ originate from the strong S77(3)i and 2/ (l)i dynamics, respectively. The 
£7(l)i force is attractive in the <it> channel but repulsive in the <bb> channel, such that 
the top but not the bottom acquires dynamical mass under the condition 



Vb < Vcxit 



/8tt 2 



< yt 



(17) 



Equivalently, this implies that the composite Higgs $ t but not develops a VEV, i.e., 
Vt 7^ and v b = 0, in contrast to the simplest 2HDE of top-condensate model analyzed 
in the previous subsection (where v tib 7^ 0). For /x < A, the composite Higgs doublets <3> t 



111 unvj ir Kst~i\^\j u J. v/ii y v v nvi W Z O / J ' / 

and develop the gauge invariant kinetic terms 
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= Z„ t {D^ t )\D^ t ) + Z, 6 (iy$ 6 )t(D M $ 6 ) , with Z* tb = §^V 2 t, b ^ , (18) 

i i 

as well potential term V(Z£ <J>^, $&). 

From ([iTf), we note that unless «i is unnaturally large (compared to «), both j/t and 
Ub should be close to the critical value y crit = y / 87r 2 /3 ~ 5.13 at the scale A. At the 
lower energy scale /x(< A), y b is still close to y t , based on the RG analysis shown in 
Fig. [7[ From the RG evolution, we find t/[,(m t ) = 2.7 — 2.1 and yt{m t ) = 2.9 — 2.2 for 
A = 1 — 10 TeV, with the typical boundary conditions 2/t(A) = 5.5 and y&(A) = 4.5 . The 
precise boundary values of yt,b(A) depends on the detailed dynamics of topcolor breaking 
via the parameters («, But for (/«, «i) not much larger than 0(1), 2/&,t(A) should be 
reasonably close to the critical value y cr i t at the scale A. Furthermore, as shown in Fig. |7], 
the infrared behavior of yt,&(/x) a t the weak scale is not sensitive to the possible variations 
of their boundary values at the scale A. Therefore, this model generically predicts a large 
yb of 0(2 — 3) at the weak scale. 



Another essential feature of the TCATC model |T0[ is that the topcolor interaction 
must not be responsible for the whole EWSB, but is mainly responsible for the top quark 
mass generation. As a result, the dynamical scale can be as low as A = 0(1) TeV (which 
avoids the severe fine-tuning needed in the minimal models pji~2"l) and correspondingly, 



vt — 64 — 97 GeV (for A = 1 — 10 TeV) according to the Pagels-Stokar formula [[3^ 




< = V^ m i M In ^7 + c (19) 



where cq = 0(1) is a constant. The b quark gets large portion of its mass from topcolor 
instanton effects ||10|| . The EWSB is mainly driven by the usual extended technicolor 
(ETC) [[37],[38|] (or the equivalent Higgs) interaction which gives small masses [~ O(GeV)] 
to all fermions (including t and b). In addition, this model predicts three physical top- 
pions with masses around 0(m t ). The smaller vacuum expectation value v t [estimated 
by (ID in TCATC models], as compared to the full VEV (v ~ 246 GeV), makes the 
Yukawa coupling of the top to $f stronger than that in the SM, which is consistent with 
the predictions from the RG analysis in Fig. [7|. The large-Af. calculation suggests 



that the neutral components (hb, A&) of are degenerate and can be the lightest Higgs 
bosons with masses of 0(100) GeV to a few hundred GeV. Hence, this model also predicts 
light scalars which couple to the bottom quark strongly via Yukawa interaction. 

To illustrate how Tevatron (Run II) and LHC can test this model via the reac- 
tion pp/pp — > bbh{-^> bb) + X , we compare yb{n) (predicted in Fig. |7|) with the 
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model-independent bound on K min y^ derived in Sec. II. Or, equivalently, we compare 
y b {fi)/yl M {fi) with K min y^/yl M {fi) [= K min m p b ° l /m b (/i)]. The result is shown in Fig. |b. 
It is clear that the Tevatron Run II can already provide useful information on this type of 
models. The LHC data can further extend the coverage of the mass range up to 1 TeV. In 
Fig. || we present the exclusion curve at the 95% CL. and the discovery reach at the 5a 
level for the LHC with 100 fb _1 of luminosity. If we relax the mass degeneracy condition 
(suggested by the large- N c analysis) and assume that the degeneracy of h b and A b does 
not hold well enough to be within the mass resolution of the detector (cf. Sec. II), then 
the contours in Figs. [| and |9] will move up by an overall factor of v^2 for most of the mass 
region. Fig. |9] shows that even in this non-degeneracy case, the LHC (with a 100 fb _1 
luminosity) can discover the Higgs boson h b or A b with a mass up to 1 TeV at the 5a 
level, since the theory curves always lie above y b /yl M = 100 for M hb{Ab) > 50 GeV (cf. 
Figs. |b). 



C. In Dynamical Left-Right Symmetric Extension of the Top-condensate Model 



Finally, we consider the left-right symmetric extension JTTJ] of the top-condensate 
model, which postulates the gauge structure (?lr = SU(3) c ®SU(2) l ®SU{2) r ®U{1)b-l 
at a high energy scale A. This model has many attractive features. For example, parity 
violation can appear naturally via the spontaneous symmetry breaking and the known 
quarks and leptons fit economically into fundamental representations of the gauge group. 
A dynamical see-saw mechanism can also be realized in this scenario, which naturally 
yields the small neutrino masses. 

At the compositeness scale A, a set of NJL-type four- Fermi interactions are gener- 
ated, which produce a composite Higgs sector at the lower scale n(< A). The sym- 
metry breaking pattern occurs via two steps: first, Qlr breaks down to Q SM = 
SU(3) C ® SU{2)l eg) U{1)y at a scale /i = Ar; second, the remaining standard model 
group (?sm is broken down to SU(3) C ® U(l) em at the electroweak scale of O(100 GeV). 
The composite Higgs sector of this model contains a scalar bi-doublet two scalar dou- 
blets xl and xr, and a singlet scalar a, with the quantum number assignments (1, 2, 2, 0), 
(1, 2, 1, —1), (1, 1, 2, —1) and (1, 1, 1, 0), respectively. They are defined as: 

| { Xl + vl\ ( Xr + vr 

<I> - V2 .n " XL = 



V 



bl + V2 

V2 J 



V 



V2 



XR 



V 



V2 
Xr 



a 



(20) 



where the vacuum expectation value vr is much larger than the other VEVs {v i 5 2 and vl) 
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and is responsible for the first step breaking of the left-right symmetry. The true VEV 
(v ~ 246 GeV) for the EWSB is determined by 



+ vt) + \ \{v\ + vl) - - vlf + {A Vl v 2 f 



V 



2 i 2 , 2 

-V X +v 2 +v Ll 



(21) 



where the approximate relation holds, because vr ^> v x<2 , vl- Note that a nonzero vl 



(which may be relatively small) implies v 12 = yv x + v 2 < v ~ 246 GeV. Furthermore, 
the singlet scalar a does not develop VEV. 

The relevant Yukawa interactions and mass terms for the top-bottom sector can be 



written as ITT 



£f F = m t tt + m b bb + ^t{y X (j>\ + y 2 0°) t + -±b ( yi $ + y 2 </>°) b 



+ Y>l [yi<Pi - 1/202 )t R + t L [yx(f)J - y 2 <Pi) b R + h.c. 
m t = (yiV! + y 2 v 2 ) /V% , m b = {y x v 2 + y 2 v x ) / V2 



:22) 



which only involve the scalars in the bi-doublet $. To give the correct top mass, tan/3 
(= v x /v 2 ) is constrained to be in the range of 1.3 — 4.0. The formation of dynamical 
condensates or the VEVs of the composite Higgs scalars requires the Yukawa couplings 
2/i,2 to be above their critical value at the compositeness scale. Consequently, y Xj2 at the 
weak scale can be naturally large [of ~ 0(1)]. The RG analysis |ljj indeed shows that for 
A to be in the range of 10 5 to 10 19 GeV, the Yukawa coupling yi(fi) varies from about 2.1 
to 1.2 at the scale /i = 0(100 — 1000) GeV. Since the Yukawa coupling y 2 (fi) satisfies the 
same RG equation as that of yi(fi) (after interchanging y x and y 2 ) |i~Tl , the infrared value 
of y 2 (fi) is also naturally large [of 0(1)], and is not sensitive to the boundary condition 
at the compositeness scale.0 Furthermore, the mass of the neutral OP-even and OP-odd 
scalars Re(0°) an d lm(0°) m &y be as light as about O(100) GeV fH|| . We thus expect 
that measuring the production rate of these light scalar bosons via the <pbb mode at the 
Tevatron and the LHC can effectively test this model. 

Before concluding this section, we note that in the three types of dynamical models 
discussed above, the relevant composite Higgs scalars (having large Yukawa coupling yb) 
do not couple to t~t + mode at the tree level. This is in contrast to the case of MSSM 
where the lepton-Higgs Yukawa couplings (such as A-t + -t~ etc) are enhanced in the 



8 Here, the compositeness scale A can be as low as 100 TeV and the left-right breaking scale 
A R = v R is around of O(10 TeV) flj. 
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same way as that for the down-type quarks in the large tan/3 region. Therefore, further 
combining the bb(fi(—> r~r + ) channel into our analysis would be useful to discriminate the 
above dynamical models from the MSSMQ, should a signal be observed. 



IV. CONSTRAINTS ON MSSM PARAMETERS AND IMPLICATIONS FOR 
MODELS OF SOFT-BREAKING OF SUSY 



Supersymmetry (SUSY) is one of the most natural extensions of the SM, mainly be- 
cause of its ability to solve the hierarchy problem, as well as for its capacity to imitate 
the current experimental success of the SM, despite the plethora of the introduced new 
particles and free parameters [0. The Minimal Supersymmetric SM (MSSM) || requires 
a two Higgs doublet extension of the SM |J together with the corresponding supersym- 
metric partners. The model includes all renormalizable interactions that respect the 
standard gauge group SU (3)c <S> SU (2) l <8> U(l)y and supersymmetry. In order to prevent 
potentially dangerous baryon and lepton number violating interactions, invariance under 
a discrete -R-parityQ is also required. To be compatible with data, supersymmetry has 
to be broken. The breaking of SUSY is parametrized by the general set of soft-breaking 
terms, which, in principle, should be deduced from a specific underlying model for SUSY 
breaking, such as the supergravity and gauge-mediated |14| models. In this section, 
we discuss the potential of the Tevatron and the LHC to test MSSM via measuring the 
production rate of <pbb mode. We shall also discuss the implication of this result on various 
supergravity and gauge-mediated models of soft SUSY-breaking. 



A. Bottom Yukawa Couplings and the MSSM Higgs Sector 

In the MSSM, the Higgs couplings to the SM fermions and gauge bosons involve two 
new free parameters at the tree-level, which are the vacuum angle f3{= axctanv u /vd) and 
the Higgs mixing angle a. These couplings are shown in Table [V]. We see that the MSSM 



9 The r Yukawa couplings to other possible composite scalars are not yet well specified in 
the top-condensate/topcolor models, while for the dynamical left-right model the r Yukawa 
couplings are expected to be naturally small, not much different from the SM value pl| , |40| . 

10 The -R-parity is defined in such a way that SM particles are even under R and their super- 
partners are odd. 
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Higgs couplings to the gauge boson pairs are always suppressed relative to that of the 
SM, while their couplings to the down (up) -type fermions are enhanced in the large (small) 
tan/3 region. These enhanced Yukawa couplings are of great phenomenological importance 
for the Higgs detection and especially for probing the associated new dynamics in the top- 
bottom sector. Alternatively, we can choose tan/3 and the pseudoscalar mass rriA as two 
free parameters. Then, at the one-loop level, a can be calculated from 

tan 2a = tan 2/3 (m A + m| 



m A — m 7 + 



cos 2/3 



(23) 



with a G (— 7r/2,0). Here the parameter e t represents the dominant top and stop loop 
corrections which depend on the fourth power of the top mass rrit and the logarithm of 
the stop mass M~: 

3GW 4 ( Mf \ 
£t= rK o . \ J og Mr + 1 • (24) 



27r 2 sin (5 



Note that for large tan/3, the bottom and sbottom loop corrections can also be important. 
Hence, in our numerical analysis below (cf. Sec. IVB), we have included the complete 
radiative corrections with full mixing in the stop and sbottom sectors, and the renormal- 
ization group improvements are also adopted.^ As shown in Table |V|, the A-b-b coupling 
has no explicit a dependence. The bottom Yukawa couplings y^ and y^u are a- and 



/3-dependent, their magnitudes relative to the SM prediction are displayed in Fig. [TT3| a as 
a function of for various tan/3 values. It shows that for uia above ~ 110 GeV, the 
h-b-b coupling quickly decreases, approaching to the SM value for all tan/3, while the H-b-b 
coupling increases for large tan/3. Therefore, we expect that in the large tan/3 region, the 
production rate of Abb or hbb can be large for small m^, while the rate of Abb or Hbb 
are enhanced for large rriA- Whether the signals of the two Higgs scalars (A and h or A 
and H) can be experimentally resolved as two separate signals (e.g., two bumps in the bb 
invariant mass distribution) depends on their mass degeneracy. 

The MSSM Higgs boson mass spectrum can be determined by taking the second 
derivative on Higgs effective potential with respect to the Higgs fields. At tree-level, the 
resulting Higgs masses obey the relations: rrih < cos 2/3, mz < run, rrih < rriA < 
rriH, and my/ < m#±. However, these relations are substantially modified by radiative 



We have used the HDECAY program [35 to compute the Higgs masses, couplings and decay 
branching ratios. 



19 



corrections [41]. Including the important contributions from top and stop loops, the 
masses of h and H can be written as: 



m h,H = \ | ( m2 + e *) T ( m2 + e t) 2 - 4e * ( m l c °s 2 /3 + ™\ sin2 Z 3 ) - 4m 2 z m 2 A cos 2 2/3 

(25) 

where M 2 = m\ + m 2 ^ and the parameter e t is defined in (fM|). To analyze the Higgs 
mass degeneracies, we plot the mass differences m^-m^ and m^-m^ in Fig. pTj[b using the 



complete radiative corrections to the Higgs mass spectrum ||35|| . We see that for the large 
tan/3 values, the pseudoscalar A is about degenerate in mass with the lighter neutral scalar 
h below ~ 120 GeV and with the heavier neutral H above ~ 120 GeV. This degeneracy 
indicates that the <pbb signal from the MSSM generally contains two mass-degenerate 
scalars with similar couplings, and thus results in a stronger bound on tan/3 by about a 
factor of v2- 

Finally, we note that in Fig. |HJ all soft-breaking mass parameters were chosen to 
be 500 GeV. Various choices of SUSY soft-breaking parameters typically affect these 
quantities by about 10%-30%. To illustrate these effects, we plot in Fig. |TT] the same 
quantities, but changing the right-handed stop mass to Mi = 200 GeV in (a-b), and in 
(c-d) we use the "LEP2 II Scan A2" set of SUSY parameters^ for comparison. 

Because the MSSM predicts a large bottom quark Yukawa coupling for large tan/3, 
and the mass of the lightest neutral scalar has to be less than ~ 130 GeV, we expect that 
the Tevatron and the LHC can test this model via measuring the (f)bb production rate. 
In the following, we shall discuss the range of the m^-tan/3 plane that can be explored 
at various colliders. Some models of SUSY soft-breaking predict a large tan/3 with light 
Higgs scalar (s), and thus predict a large <fibb rate. Without observing such a signal, one 
can put a stringent constrain on the model. This will also be discussed below. 



12 The parameters uiq and Mi are fixed at 1 TeV, fi is chosen to be —100 GeV and rat = 
175 GeV. The scalar trilinear coupling Ai is fixed at \/6 TeV, corresponding to the maximal 
left- and right-handed top-squark mixing. Detailed prescription about this set of parameter scan 
can be found in Ref. HH- 
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B. Constraints on MSSM from ipbb production at Tevatron and LHC 



To use the model-independent result of K min obtained in Sec. II to constrain the m^- 
tan/3 plane in the MSSM, one needs to calculate the SUSY Higgs boson masses, decay 
branching ratios, and their couplings to the bottom quark for a given set of the soft 
breaking parameters. In the following numerical analysis, we use the HDECAY code to 
include the full mixings in the stop/sbottom sector with QCD and electroweak radiative 
corrections fl35| . For simplicity, we assume that the superpartners are all heavy enough 
so that the decays of the Higgs bosons into them are forbidden. Under this assumption, 
we find that the decay branching ratio of h — > 66 is close to one for the relevant region of 
the parameter space. 

As explained above, we combine signals from more than one scalar boson provided 
their masses differ by less than Am exp , which is the maximum of the experimental mass 
resolution (cf. footnote-7) and the natural decay width of the scalar boson. Since the 
results of Sec. II are given in terms of the minimal enhancement factor K m[n defined in 
(1), we need to convert them into exclusion bounds in the m^-tan/5 plane of the MSSM, in 
case that a signal is not found. The bound on tan fi min (with the possible mass degeneracy 
included) can be derived from that on K min (for a single scalar) by requiring 



=h,H 




tan 2 /? BR{A (A™ex P - \AM A<t> \) BR(0 - 66) > K 2 min , (26) 



where y^ M and yf denote the b quark Yukawa coupling in the SM and the MSSM (with 
<f> — h or H), respectively. Inside the argument of the ^-function, AM^ is the mass 
difference between A and 0. Thus, the equality sign in the above relation determines the 
minimal value tan/3 m i n for each given if mm- 

To estimate the exclusion regions in the m^-tan/5 plane, a set of soft breaking pa- 
rameters has to be chosen, which should be compatible with the current data from the 
LEP II and the Tevatron experiments, while not much larger than 1 TeV. For simplicity, 
we choose all the soft SUSY breaking parameters (and the Higgs mixing parameter-/^) 



to be 500 GeV as our "default" values, i.e., M so ft = 500 GeV. In Fig. |12|a, we show 
the 95% CL. exclusion contours in the m^-tan/3 plane derived from the measurement of 
pp/pp — > 066 — ► 6666, using this "default" set of SUSY parameters. The areas above the 
four boundaries are excluded for the Tevatron Run II with the indicated luminosities, and 
for the LHC with an integrated luminosity of 100 fb _1 . Needless to say, different choice 
of SUSY parameters, such as the mass and the mixing of the top squarks and the value 
(and sign) of the parameter /i, would modify this result. To compare the potential of the 
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Tevatron and the LHC in constraining the MSSM parameters via the <pbb {(f) = h, A, H) 
production to that of the LEP II experiments via Z<p and hA production, we consider 
one of the "benchmark" parameter scans discussed in p5| , |4*3|| , which is called the "LEP II 
Scan A2" 12 set. For this set of SUSY parameters, the LEP II exclusion contour 0,^] at 
the 95% CL. is displayed in Fig. [I2|b, for a center-of-mass energy of 200 GeV and an inte- 



grated luminosity of 100 pb -1 per LEP II experiment. As shown in Fig. |l^b, the Tevatron 
Run II result, in comparison with the LEP II result, can already cover a substantial region 
of the parameter space with only a 2 fb _1 luminosity. Thus, detecting the cj)hb signal at 
hadron colliders can effectively probe the MSSM Higgs sector, especially for models with 
large tan (3 values. Furthermore, for uia ~ 100 GeV, Tevatron Run II is complementary 
with LEP II, because the latter is not sensitive to that region of parameter space. The 
LHC can further probe the MSSM down to tan/3 ~ 7 (15) for m A < 400 (1000) GeV. 
This is also shown in Fig. O using the "default" set of SUSY parameters, in which the 
region above the upper curve is the discovery contour at the 5a level for the LHC with 
an integrated luminosity of 100 fb -1 , and the area above the lower curve can be excluded 
at 95% C.L, if a signal is not found. 

For completeness, we also present the exclusion contours in the ra^-tan/? and 
planes for both the "default" and the "LEP II Scan A2" sets of SUSY parameters. They 
are shown in Figs. and [15|. Again, we see that the Tevatron Run II and the LHC bounds 
sensitively and complementarily cover the MSSM parameter space in contrast with the 
LEP II results. We have also studied the bounds with the "LEP II Scan Al" inputs |42]j4^ 
and found that the exclusion contours from the <pbb production are similar to those with 
"LEP II Scan A2" inputs. The most noticeable difference is that the theoretically allowed 
range for rrih becomes smaller by about 10 GeV in the "Scan Al" set, as compared to the 
"Scan A2" inputs. 

Even though our analyses, described above and in Sec. II, are quite different from that 
of Ref. ||, the final bounds at the LHC happen to be in qualitative agreement. We also 
note that our bounds on the m^-tan/? plane improve considerably the one obtained in 
Ref. ||44|| , in which the pp — > <pbb — > r + r~bb production rate at the Tevatron Run I data 
was compared to the MSSM prediction. Though, we do not choose to explicitly present 
projected results for the Tevatron Run I data, we encourage our experimental colleagues 
to pursue this analysis on the existing Run I data sample, as it seems likely that one could 
obtain useful information even with the lower luminosity and collider energy of Run I as 
well somewhat lower 6-tagging efficiency 

Before concluding this subsection, we remark upon the effects on our bounds from 
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the possible radiative corrections to the <j)bb production process. As mentioned earlier 
in Sec. II, one of the dominant correction is from the next-to-leading order (NLO) QCD 
loops, which are not currently available for the (ftbb signal and background cross sections. 
However, aside from the QCD corrections to the (ftbb vertex (part of that can be simply 
included into the running of the <j)bb Yukawa coupling or the running 6-mass), there 
are pentagon loops formed by the virtual gluons radiated from an initial state quark 
(gluon) and re-absorbed by the final state b quark with the <pbb vertex included in the 
loop. Such QCD corrections are not factorizable so that a consistent improvement of our 
results is impossible before a full NLO QCD analysis is completed^. Putting aside the 
complexity of the full NLO QCD corrections, we briefly comment upon how the radiative 
corrections to the running (fibb Yukawa coupling affect our final bounds. The well-known 
QCD correction, Eq. fljjD , alone will reduce the running mass m^/i) by about 40% from 
the scale /i = ml° le ~ 5 GeV up to the weak scale of 0(200) GeV (cf. the solid curve of 



Fig. |16|). This is however not the full story. The complexity comes from the finite SUSY 



threshold correction in the large tan/? region which are potentially large [47-49] . In this 



case, as shown in Ref. |[48|| , the dominant one- loop SUSY correction contributes to running 
6-mass a finite term so that m^fi) at the SUSY scale /i = /zr = M so f t is multiplied by 
a constant factor 1/[1 + A5(SUSY)] , which appears as a common factor in the bottom 
Yukawa couplings of all three neutral Higgs bosons. For large tan/?, Ab(SUSY) contains 
the following tan/3-enhanced terms from sbottom-gluino and stop-chargino loopsQ, 

At(SUSY) = + 

V fUh J \ rrih J 



-/i tan/3 



16vr 2 13 
with the function F defined as: 



(27) 

glm~ g F (m b , m~ b , m- g ) + y t A t F (m h m h \i) - g 2 2 M 2 F {m h m 2 , /i) 



xy In x/y + yz In y/z + zx In zjx 
(x - y)(y - z)(z - x) 



13 



11 



This point has also been recently discussed in Ref. 
Such a full NLO QCD calculation is beyond the scope of our current study. A systematic 



calculation for this is in progress [46]. 



15 We thank K. Matchev for discussing his published results in Ref. [48|, and to him and 
W.A. Bardeen for discussing the issue of the running 6-mass. Our convention of the MSSM 
Higgs parameter differs from that of Ref. [E3] by a minus sign. 
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In these equations, the MSSM Higgs parameter /i should not be confused with the usual 
renormalization scale //r. In (pTD , we have assumed, for simplicity, mass degeneracy for 
the top and bottom squarks, respectively. Eq. (|27| ) shows that the SUSY correction to 
the running rrib is proportional to tan/3 and fi. Thus, this correction is enhanced for large 
tan/3 and non-negligible in comparison with the QCD corrections. Also changing the sign 
of fi will vary the sign of the whole correction Ab(SUSY) and implies that the SUSY 
correction can either increase or decrease the running 6-mass at the energy scale around 
of 0(M so f t ). Normally, when defining the running coupling using the Collins- Wilczek- 
Zee (CWZ) scheme fBTlfl , only the /i^-dependent contributions are included while all the 



/^-independent terms are absorbed into the corresponding Wilson coefficient functions. 
However, since the /i^-independent contribution A;, (SUSY) is not small for large tan/? 
and a full NLO SUSY calculation is not yet available, we include A;, (SUSY) to define an 
"effective" running coupling/mass of 6-quark even below the SUSY threshold scale M so f t . 
This could give a rough estimate on the large SUSY loop corrections from the <fi-b-b vertex. 
Obviously, when fiR is much smaller than M so f t , the CWZ scheme should be used. Hence, 
in Fig. [16|, we only show the "effective" running mass m;,(/i) down to about 100 GeV which 



is the relevant energy scale and the mass scale (~ Mmggs) considered in this paper. Fig. [1C 
illustrates that due to the SUSY correction the "effective" running b mass at the weak scale 
can be either larger or smaller than the SM QCD running value (~ 3 GeV), depending 
on the choice of the sign of /i parameter (and also other soft-breaking parameters such as 
the trilinear coupling A t and masses of the gluino, gaugino, stop and sbottom). For the 
"default" set of SUSY parameters used in our analysis, all the soft-breaking parameters 
are set to 500 GeV for simplicity. It happens to be the case that the SM QCD and SUSY 
corrections nearly cancel each other so that the "effective" running mass rrib is very close to 



the pole mass value (~ 5 GeV) for the scale above ~ 100 GeV (cf. upper curve of Fig. [16]). 
For comparison, in our analysis using the "LEP II Scan A2" set of SUSY parameters, the 
SM QCD and SUSY corrections do not cancel and tend to reduce the "effective" running 
6-mass or the Yukawa coupling ?/b(/i) which results in a weaker bound for the Tevatron 
Run II and the LHC, as shown in Fig. [l2|b.[^] The difference in the exclusion contours 
shown in Figs. |T2| a and |T2|b, derived from the measurement of the <pbb production rate 
at hadron colliders, is mainly due to the difference in the "effective" running coupling or 
mass (including the QCD and SUSY corrections), as described above.0 We therefore 
conclude that a full NLO QCD calculation is important for a consistent improvement of 



16 A detailed analysis of these effects at the Tevatron Run II is currently underway [45 



1 From Fig. [T^, it is also clear that, for the LHC bounds, the SUSY correction has much less 
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our current analysis. 



C. Interpretation of results for Models of soft breaking parameters 

The MSSM allows for a very general set of soft SUSY-breaking terms and thus is 
specified by a large number of free parameters (~ 124 though only a complicated 
subset of this parameter space is consistent with all current experimental results. It 
is therefore important to understand the mechanism of supersymmetry breaking (which 
presumably occurs at a high energy scale joTR ) and to predict the soft parameters at the 
weak scale from an underlying model. Many alternative ideas about how supersymmetry 
might be broken, and how this will result in the low energy soft breaking parameters 
exist in the literature, including the supergravity inspired (SUGRA) models and gauge- 
mediated SUSY breaking (GMSB) models. In this section we examine the sensitivity of 
the (pbb process to probe a few models of SUSY breaking, concentrating for the most part 
on the popular SUGRA and GMSB models. However, there are also other interesting 
ideas to which the <ftbb process may provide interesting information, because these models 
naturally prefer a large tan/5. A few examples include the SO (10) grand unification 
theories [0] (SO(10) GUTs); the infrared fixed-point scenario [|53|; and also a scenario 



with compositeness, the "more minimal supersymmetric SM" |54 . 



1. Supergravity Models with large tan/3 



The supergravity inspired (SUGRA) models |55| incorporate gravity in a natural man- 
ner, and solve the problem of SUSY breaking through the introduction of a hidden sector, 
which breaks SUSY at a very high scale [~ O(10 n ) GeV] and interacts with the MSSM 
fields only gravitationally. This model offers an exciting glimpse into the possible connec- 
tion between the heavy top quark and the EWSB by allowing radiative breaking of the 
electroweak symmetry, in which the large top quark Yukawa coupling can drive one of the 
Higgs masses negative at energies ~ rriz- In the limit of large tan/3 the bottom and tau 



Yukawa couplings can also play an important role [5~o] . 

Under the assumption that the gravitational interactions are flavor-blind, this model 
determines the entire SUSY spectrum in terms of five free parameters (at the high energy 



impact since the relevant tan/3 values become much lower, around of 0(2 — 15). 
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scale of the SUSY breaking) including a common scalar mass (too), a common gaugino 
mass (M), a common scalar tri-linear interaction term (A), tan (3, and the sign of the Higgs 
mixing parameter (sgn (//)). The weak scale particle spectrum can then be determined by 
using the renormalization group analysis to run the sparticle masses from the high scale 
to the weak scale. 

Though large tan/5 is not required by the minimal SUGRA model, it can naturally 
be accommodated, as demonstrated in > "where it was found that a large tan (3 also 

generally requires that the pseudoscalar Higgs mass be light (tua < 200 GeV for tan (3 > 
30), because the enhanced b and r Yukawa couplings act through the renormalization 
group equations to reduce the down-type Higgs mass term at the weak scale, thus resulting 
in a light Higgs spectrum. Since the importance of the b and r effects in the RG analysis 
increases with larger tan /3, as tan j3 increases the resulting tua decreases, making the large 
tan p scenario in the SUGRA model particularly easy to probe through the (pbb process. 
From the limits on the m^-tan /3 plane derived above in Section IVB, it thus seems likely 
that from the data of the Tevatron Run II with 2 fb _1 of integrated luminosity, a large 
portion of the minimal SUGRA model with tan (3 > 20 may be excluded. 



2. Gauge-mediated SUSY Breaking Models with Large tan/3 



Models with GMSB break SUSY at a scale which is typically much lower than that 
present in the SUGRA models. The supersymmetry is generally broken in a hidden sector 
which directly couples to a set of messenger chiral superfields. This induces a difference 
in mass between the fermion and scalar components of the messenger fields, which in turn 
generates masses for the gauginos and sfermions of the MSSM fields via loops involving 



the ordinary gauge interactions A generic feature of this scenario is that because 

of the relatively low scale of SUSY breaking, the gravitino acquires a much smaller mass 
than in the SUGRA scenarios, and is generally the lightest super symmetric partner (LSP). 
While specific models of GMSB vary as to their assumptions and relevant parameters, 
generally what must be assumed is the field content of the messenger sector (including 
transformation properties under the gauge group and number of multiplets in the theory) 
and the scale at which SUSY is broken in the hidden sector. 

The minimal GMSB models can also result in a radiative breaking of the EWSB, 
through the renormalization group evolution of the Higgs masses (driven by the large 
top Yukawa coupling) from the effective SUSY breaking scale to the weak scale. As in 
the SUGRA model case, this evolution can drive the mass term of the up-type Higgs 
negative at the weak scale, thus breaking the electroweak symmetry. In fact, because 
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the effective SUSY breaking scale in a GMSB model is typically much lower than in the 
SUGRA model (and thus closer to the weak scale), in order for the proper EWSB to 
occur, it was demonstrated in |50] that a large tan/3 is required (about 30-40). However, 
because the effective SUSY breaking scale is typically much lower than in the SUGRA 
model, the large effects of the b and r on the Higgs mass running do not reduce the 
Higgs spectrum to the degree that occurs in the SUGRA model with large tan/3, and 
thus result in a heavier A with mass of about 400 GeV. So, this particular model would 
only be explored through the <pbb process at the LHC. However, more general analyses of 
the GMSB scenario |6l-|S^] , introducing more degrees of freedom in the messenger sector 
than the minimal model, can allow for large tan/3 and relax tua to be as low as about 200 
GeV. Thus, through the <frbb production these more general models can be first probed at 
the Tevatron for the relevant mass range, and then largely explored at the LHC. 



V. CONCLUSIONS 



It remains a challenging task to determine the underlying dynamics of the electroweak 
symmetry breaking and the flavor symmetry breaking. Either fundamental or composite 
Higgs boson(s) may play a central role in the mass generation of the weak gauge bosons 
and the fermions. The heavy top quark, with a mass on the same order as the scale of 
the electroweak symmetry breaking, suggests that the top quark may play a special role 
in the mechanism of mass generation. In this work, we have shown that in the typical 
models of this type, the bottom quark, as the weak isospin partner of the top quark, can 
also participate in the dynamics of mass generation, and serves as an effective probe of 
the possible new physics associated with the Higgs and top sectors. 

We have presented a model-independent analysis on Higgs boson production in as- 
sociation with bottom quarks, via the reaction pp/pp — > <pbb — > bbbb, at the Tevatron 
Run II and the LHC. We have computed the QCD bbbb and bbjj backgrounds, and the 
electroweak Zbb background to illustrate that it is possible to extract the signal from 
these large backgrounds by employing a suitable search strategy. The scale and the PDF 
dependencies in our signal and background calculations are also examined and they in- 
dicate that the NLO QCD corrections to the production rate could be large, and thus 
including the complete NLO QCD corrections in future improvements will be very useful. 
Using the complete tree level calculation with an estimated QCD fc-factor of 2, we derive 
the exclusion contour for the enhancement factor (in the coupling of <pbb relative to that 
of the SM) versus the Higgs mass at the 95% C.L., assuming a signal is not found. 

We apply these results to analyze the constraints on the parameter space of both the 
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composite models and the MSSM (in the large tan/3 region) with naturally large bottom 
Yukawa couplings. For the composite Higgs scenario, we first consider the two-Higgs- 
doublet extension of top-condensate model and then analyze the topcolor model, where 
the b quark Yukawa couplings are naturally large due to the infrared quasi-fixed-point 
structure and the particular boundary conditions for (y b , y t ) at the compositeness scale. 
Our analysis shows that the Tevatron Run II with a 2 fb _1 of luminosity can exclude 
the entire parameter space of the simplest two-Higgs-doublet extension of top-condensate 
model, if a signal is not found. For the topcolor model, the Tevatron Run II is able to 
detect the composite Higgs h b or A b up to ~ 400 GeV and the LHC can extend the mass 
range up to ~ 1 TeV. Similarly, this production mechanism can be used to effectively test 
the dynamical left-right symmetric model. 

To confirm the MSSM, it is necessary to detect all the predicted neutral Higgs bosons 
h, H, A and the charged scalars H ± . From LEP II, depending on the choice of the MSSM 
soft-breaking parameters, the current 95% C.L. bounds on the masses of the MSSM Higgs 



bosons are about 70 GeV for both the CP-even scalar h and the CP-odd scalar A 651. It 



can be improved at LEP II with higher luminosity and maximal energy, but the bounds 
on the Higgs masses will not be much larger than ~ mz for an arbitrary tan/? value. The 
Wh and WH associated production at the Tevatron Run II can further improve these 
bounds, if a signal is not observed. At the LHC, a large portion of parameter space can 
be tested via pp — ► tt+h(— > 77) +X, and pp — > h(—> ZZ*) +X, etc p6| , |67| . A future high 
energy e + e~ collider will fully test the MSSM Higgs sector [55E5| through the reactions 
e + e _ — > Z + h(H), A + h(H), H + H~ , etc. In this paper, we demonstrate that studying 
the (f)bb channel at hadron colliders can further improve our knowledge on MSSM. The 
exclusion contours on the m^-tan/? plane of the MSSM shows that the Tevatron and the 
LHC are sensitive to a large portion of the parameter space via this mode. It therefore 
provides a complementary probe of the MSSM Higgs sector in comparison with that from 
LEP II. The implications of these bounds in the parameter space on both the supergravity 
and the gauge-mediated SUSY breaking models are further discussed. We find that the 
<f)bb process can effectively test the models with either scheme of the SUSY soft-breaking 
in the large tan/3 scenario. 

In conclusion, the stringent constraints obtained by studying the associated <ftbb pro- 
duction at hadron colliders for both the composite and the supersymmetric models show 
the utility of this search mode. However, much work remains to be done. For exam- 
ple, a 6-trigger would be essential for this analysis. The CDF group at the Fermilab has 
demonstrated that it is possible to detect events with four or more jets and two or more b- 
tags [70], which can be significantly improved with the implementation of 6-trigger at the 
Run II of the Tevatron |7T|. However, the large QCD 4-jet background rate at the LHC 
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can potentially make triggering on the bb<f>(— > bb) events difficult, though it is expected 
that a 6-trigger would become more efficient for a heavier (pseudo-) scalar 0. (This is be- 
cause the 6-jets from the decay of a heavier are more energetic, and its QCD background 
rate drops rapidly as a function of the transverse momentum of the triggered jet.) We 
hope that the interesting results afforded by studying this channel will stimulate interest 
of our experimental colleagues in working on these problems. 

We have found that the (j)bb process complements Higgs searches in other channels, 
and thus it is expected that experimental searches for this signature at the Tevatron 
Run II (and possibly beyond) and the CERN LHC will provide interesting and important 
information about the mechanism of the electroweak symmetry breaking and the fermion 
mass generation. 
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TABLES 



Process 


Acceptance Cuts 


Pt Cuts 


AR Cut 


AM Cut 


<f>bb 


4923 


1936 


1389 


1389 


Zbb 


1432 


580 


357 


357 


bbbb 


5.1 x 10 4 


3760 


1368 


1284 


bbjj 


1.2 x 10 7 


1.5 x 10 6 


6.3 x 10 5 


5.9 x 10 5 



TABLE I. The signal and background events for 2 fb of Tevatron data, assuming 
77i0 = 100 GeV, 2A 7770 = 26 GeV, and K = 40 after imposing the acceptance cuts, px cuts, 
and reconstructed cuts described in the text. (A fc-factor of 2 is included in both the signal 
and the background rates.) 



7770 (GeV) 


p^ ] Cut 


P { P Cut 


75 


35 


25 


100 


50 


30 


125 


60 


35 


150 


70 


45 


175 


85 


55 


200 


90 


60 


250 


125 


80 


300 


150 


105 


350 


175 


190 


400 


190 


120 


500 


245 


160 


800 


390 


260 


1000 


500 


320 



TABLE II. The optimal p T and p T cuts for isolating a Higgs boson of mass from the 
QCD bbbb background. 
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Process 


2 or more 6-tags 


3 or more 6-tags 


4 6-tags 


4>bb 


1139 


660 


180 


Zbb 


293 


170 


46 


bbbb 


1054 


610 


166 


bbjj 


1.2 x 10 5 


2141 


4 


Significance 


3.3 


12.21 


12.25 



TABLE III. The signal and background events for 2 fb _ of Tevatron data, assuming 
777,0 = 100 GeV, 2At770 = 26 GeV, and K = 40 for two or more, three or more, or four 6-tags, 
and the resulting significance of the signal. 



777 (GeV) 


N s 


Tevatron 

N B 


N s 


LHC 


N B 


75 


583 


640 


3.4 xlO 6 




4.8 xlO 6 


100 


180 


216 


2.0 xlO 6 




3.0 xlO 6 


150 


58 


92 


9.2 xlO 5 




1.2 xlO 6 


200 


17 


31 


4.2 xlO 5 




5.6 xlO 5 


250 


4.8 


8.8 


1.9 xlO 5 




2.0 xlO 5 


300 


1.3 


2.1 


83000 




70000 


500 






12000 




5700 


800 






1500 




406 


1000 






407 




70 



TABLE IV. Event numbers of signal (Ns), for one Higgs boson, and background (Nb) for a 
2 fb _1 of Tevatron data and a 100 fb _1 of LHC data, for various values of ttt,^, after applying the 
cuts described in the text, and requiring 4 6-tags. An enhancement of K = 40 is assumed for the 
signal, though the numbers may be simply scaled for any K new by multiplying by (K ncw /40) 2 . 
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Higgs 


A 


h 


H 


yu/yl u 


cot (3 


cos aj sin (3 


sin a/ sin/3 


y D /y s D M 


tan/3 


— sin a/ cos (3 


cos a / cos /3 


mv/y^v 





sm((3 — a) 


cos(/3 — a) 



TABLE V. Comparison of the neutral MSSM Higgs couplings to up-type (U = u, c, t) and 
down-type (D = d, s, b; e, fi, r) fermions and to the gauge-boson (V = W, Z) pairs. The ratios 
to the corresponding SM couplings are shown, which are determined by angles (3 and a at the 
tree-level. 
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FIGURES 




FIG. 1. Representative leading order Feynman diagrams for (f>bb production at a hadron 
collider. The decay <j> — > bb is not shown. 




FIG. 2. Representative Feynman diagrams for leading order Zbb production at a hadron 
collider. The decay Z — > bb is not shown. 
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FIG. 3. Representative leading order Feynman diagrams for QCD bbbb production at a 
hadron collider. 
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FIG. 4. Representative leading order Feynman diagrams for QCD bbjj production at a 
hadron collider. 
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O I i l-i-T ii-m-i-l 
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Pt <2> (GeV) 

FIG. 5. Figure (a) shows the distribution of the QCD bbbb background (solid curve) and 
K = 40 <pbb signal (dashed curve) cross sections in pj) at the Tevatron Run II after the ac- 

(2) 

ceptance cuts. [K = yb/(yb)sM, cf. Eq. (1).] Figure (b) presents the distribution of p T at the 
Tevatron Run II after applying the cut to , illustrating the utility of asymmetric cuts on pj) 

(2) 

and p T in extracting the <j)bb signal from the QCD background. 
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FIG. 6. (a). The model-independent minimum enhancement factor, K m - m , excluded at 95% 
C.L. as a function of scalar mass (m^) for the Tevatron Run II with 2 fb _1 (solid curve), 10 
fb" 1 (dashed curve and 30 fb" 1 (dotted curve), (b). The same factor, K m { n , excluded at 95% 
C.L. (solid curve) and discovered at 5a (dashed curve) as a function of for the LHC with 
100 fb _1 . In the above, the natural width of the scalar (1^) is assumed to be much smaller than 
the experimental mass resolution. 
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FIG. 7. Renormalization group running of y< and ?/& in the 2HDE of top-condensate model 
and the TCATC model. For the 2HDE of top-condensate model, two sets of curves are shown: 
the upper two curves (solid for yt and dotted for y&) are for the typical boundary condition 
yt(10 15 GeV) = y;,(10 15 GeV) = 6> 1 and they are too close to be distinguishable; the lower 
two curves are for the boundary condition y t (10 15 GeV) = 4, ?/{,(10 15 GeV) = 1 . In both cases, 
yt and y^ have very similar infrared values, of 0(1), at the weak scale. 



42 





I I I I I I I I I I I I I I I, 

(a) ; 


I I I I I I I I I I I I I I 

; (b) ; 






Theor. Upper Limit ! _ 
(A=10 15 GeV) I 






10 2 


' >■ 

~ tan/?=m t /m b /^^f^ 


-^A^ToTeV / / 
2fb/ / : 


10 2 


Q2_ 
C 

o 

HI — ' 


2fb>^ ^^/^^ 
/ Tevatron j 




to ^ 
>n 
\ 

>s 


10 


^/30fb' 1 : _ 

s H 


/ / Tevatron 


10 




LHC : - 








'"ioofb -1 i 






1 


III III. 





1 



100 200 200 400 

Mp(GeV) M hb (GeV) 



FIG. 8. 95% C.L. discovery reach of the Tevatron Run II and the LHC for (a) the 2HDE of 
top-condensate model and (b) the TCATC model. Regions above the curves can be discovered. 
The top curves in (b) indicate yj,(fi = m t ) values for various topcolor breaking scale A, which 
are based on the RG running analysis (cf. Fig. 7). 
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FIG. 9. The 5a discovery and 95% C.L. exclusion contours for yb{^) / yf M {n) as a function 
of M hb in the TCATC model, at the LHC with 100 fbr 1 luminosity. 
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FIG. 10. Bottom Yukawa couplings to the MSSM Higgs bosons and the mass differences, 
niA — mii and run — niA, as a function of uia for tan/3 values: 2.0, 5.0, 10.0, 20.0, 30.0, 40., 50.0. 
In (a), ybbh is m solid and ybbH is m dashed, and tan/3 decreases from top to bottom curves. 
In (b), TtiA — rrih is in solid, mu — rnA is in dashed, and tan (3 increases from the top to bottom 
curves. Here, all the SUSY soft-breaking mass parameters are chosen to be 500 GeV. 
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FIG. 11. The same as the previous figure, but in (a)-(b), we change the right-handed stop 
mass to 200 GeV, and in (c)-(d), we use the "LEP II Scan A2" set of SUSY parameters. 



46 



I i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i L 

100 200 , . 300 400 

M A (GeV) 




100 200 , . 300 400 

M A (GeV) 

FIG. 12. 95% C.L. exclusion contours in the m^-tan f3 plane of the MSSM. The areas above 
the four boundaries are excluded for the Tevatron Run II with the indicated luminosities, and 
for the LHC with an integrated luminosity of 100 fb . The soft SUSY breaking parameters 
were chosen uniformly to be 500 GeV in Fig. (a) , while the inputs of the "LEP II Scan A2" are 
used for the Fig. (b) in which LEP II excludes the left area of the solid curve. 
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FIG. 13. Discovery and exclusion contours in the m^-tan (3 plane of the MSSM for the 
LHC with an integrated luminosity of 100 fb _1 . The area above the lower boundary is excluded 
at 95% C.L., while the upper boundary is the 5a discovery contour. The soft SUSY breaking 
parameters were chosen uniformly to be 500 GeV. 
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FIG. 14. 95% C.L. exclusion contours in the m/j-tan f3 plane of the MSSM. The areas above 
the four boundaries are excluded for the Tevatron Run II with the indicated luminosities, and 
for the LHC with an integrated luminosity of 100 fb _1 . LEP II can exclude the area on the 
left-hand side of the solid curve in the lower plot. 
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FIG. 15. 95% C.L. exclusion contours in the rn^-m^ plane of the MSSM. The shaded 
areas indicate the excluded regions for the Tevatron Run II with integrated luminosities, 2, 10, 
30 fb" 1 , and for the LHC with an integrated luminosity of 100 fb _1 , as those in the previous 
figures. LEP II can exclude the area on the left-hand side of the solid curve in the lower plot. 
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FIG. 16. The running of the bottom quark mass mb(nii) as a function of the renormalization 
scale ixr. The solid curve shows the QCD evolution alone. The dashed curve further includes 
the supersymmetric corrections to the "effective" running mass, for tan/3 = 30. All soft SUSY 
breaking parameters have been fixed as 500 GeV. The dotted curve includes the SUSY corrections 
but with the sign of the Higgs-mixing parameter (i flipped. 
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